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Abstract. We use a new technique to extract the spectrum of a supernova from that of the contaminating background of its 
host galaxy, and apply it to the specific case of high-redshift Type la supernova (SN la) spectroscopy. The algorithm is based 
on a two-channel iterative technique employing the Richardson-Lucy restoration method and is implemented in the IRAF code 
specinholucy. We run the code both on simulated (SN la at z = 0.5 embedded in a bright host galaxy) and observed (SNe la at 
various phases up to z = 0.236) data taken with VLT+FORS 1 and show the advantages of using such a deconvolution technique 
in comparison with less elaborate methods. This paper is motivated by the need for optimal supernova spectroscopic data 
reduction in order to make meaningful comparisons between the low and high-redshift SN la samples. This may reveal subtle 
evolutionary and systematic effects that could depend on redshift and bias the cosmological results derived from comparisons 
of local and high-z SNe la in recent years. We describe the various aspects of the extraction in some detail as guidelines for the 
first-time user and present an optimal observing strategy for successful implementation of this method in future high-z SN la 
spectroscopic follow-up programmes. 

Key words, supernovae: general - supernovae: individual: SN 2002bo, SN 2002go, SN 2002gr - instrumentation: spectro- 
graphs - methods: data analysis - techniques: spectroscopic 



' 1. Introduction (z > 1) redshifts - the Hubble Higher-z Supernova Search 

Jzd . I Riess et all 12004) - as well as measure second order effects 



C3 The initial claim made by tw o independent team s - the (flMJ equation . of . state parame ter u = pipe 2 ) at intermediate 

High-Z Supernova Search Team fechmidt et alj l l998h and the (0,2 <,z<, 0.8) redshifts - such as ESSENCE ( MiknaitistTaTI 

Supernova Cosmology Project j Perlmutter et alj | l999|) - that 2005; Mathe s0 n et ZEjoj and the CFHT SunerNo va Legacy 

the apparent dimming of Type la supernovae (SN la) at red- s Qr SNLS (Pain & SNLS Collaboration 2QQ2) 

shifts of z ~ 0.5 implies a present accelera ting expansion of However the h sical arameters overnin the SN la ex 
the universe (R iess et alll998allPerlmutter et alJl 999). and its , . oweve ^ 1 e P ysica parameters governing t e a ex- 

— " : — — TTT — ; . . L_ 1 n plosion mechanism are today not fully understood and the em- 
subsequent confirmation with improved precision ( Tonrv et al. ... . . „ , . „"»., T , , , 

I-,,,, -I ir; — — nu„„->| VT, — • , IU™1 . — — — pineal parametnsation of the entire S N la class - the so-called 

2003; Knorj et al. 2003; Barns et al. 2004), have prompted an f™.,,. , . „, u. .... — 1|, „„„| , ■ • ,, • 

. — — , . . — r — i - — , ^ , , Phillips relation (see Philhrjs 1993) - has no physical basis, 

increased interest in these astrophysical events. On the theo 



retical side, a lot of effort and computing time has been de- 
voted to multi-dimensional modelling of the explosion in or- 



The aforementioned observational programmes will detect sev- 
eral hundreds of high-z SNe la, at which stage the statistical er- 
rors in the sample will have reached the systematic error floor. 



der to provide physical input parameters for spectral synthe „. , , , 

, , . , i TT .„ \ — r ~ „ T . ■ ■ ET ~ ~ ^ 1 ^ The cosmological effect detected via observations of Type la 

sis calculations (see Hillebrandt & Niemevei 2000 for a re- . . , , , . r — ; — ni-,, ' ' i „ 

. _ , — : : — ~ — — : — — , supernovae is indeed a subtle one (see Leibundgut 2001 for 

view). On the observational side several teams are cunently 



using SNe la both to probe the decelerating expansion at higher 



a review) and the field is limited by the systematic enors in- 
volved in such measurements and their possible evolution with 
Send offprint requests to: S. Blondin redshift. 

* Partly based on observations collected at the European Southern Spectroscopy is an ideal way to probe these potential evo- 
Observatory, Chile (ESO Programme 170.A-0519). lutionary effects through systematic comparison of high-z SN 
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la spectra with local templates. This requires the extraction 
of clean SN la spectra devoid of background contamination, 
mainly by the host galaxy and sky lines (for ground-based 
observations). For high signal-to-noise (hereafter S/N) cases, 
when the supernova is bright with respect to its immediate un- 
derlying background, standard extraction software will in gen- 
eral work well. However, when the supernova is faint compared 
to the background, e.g. at late phases, heavily superposed on its 
host galaxy or simply at a high redshift, a more elaborate tech- 
nique is required to ensure that we are not extracting non-SN 
flux. 

The purpose of this paper is to test such a method 
that has been implemented in the IRAF 1 code specinholucy 
(Lucv & Walsh 2003). It is based on a two-channel restoration 
algorithm that restores a point spread function (PSF)-like com- 
ponent in a 2D image and an underlying extended background 
separately. It is of wide astronomical use and has already been 
successfully applied to the restoration of poin t-source spec- 
tra in highly inhomogeneous backgrounds (see Lucv & W alsh 
2003), such as is often the case for high-z SNe la embedded 
in their host galaxy. The fact that the entire 2D spectrum of 
the background (with or without the inclusion of the restored 
point source spectrum) is restored in such an algorithm means 
that we have a firm hold on potential systematic errors asso- 
ciated with the restoration. This is done by simply comparing 
the residuals in the restored 2D spectrum with the statistical 
noise limit of the input 2D spectrum (the square root of the 
number of photoelectrons). Should these residuals fall below 
this limit then we can consider the restoration to be success- 
ful. This is a clear advantage over traditional spectral extraction 
routines where no secure hold on systematic errors is possible. 
Furthermore, the very nature of the algorithm presented in this 
paper implies that all non PSF-like components in the input 2D 
spectrum are automatically allocated to the background chan- 
nel. We are making no a priori assumptions on the nature of the 
contaminating extended background component and only rely 
on the point-like nature of the supernova's spatial profile. 

In section |2 we present the algorithm and its application 
to the case of supernovae embedded in their host galaxy. We 
dedicate sections[3]and|4]to the specific issues of the PSF spec- 
trum - or slit spread function (SSF) - used to restore the point 
source and of the spatial resolution kernel used to discriminate 
between the point source and the underlying extended back- 
ground, respectively, as they are of crucial importance to the 
proper functioning of specinholucy. In sections [5] and [6] we 
test the algorithm on simulated and observed data, respectively. 
Section0serves the purpose of comparing our technique to al- 
ternative spectral extraction methods, and we further discuss 
the advantages of using our approach alongside an optimal ob- 
serving strategy for its successful implementation in section[8] 

The need for comparisons between the high-z and the lo- 
cal SN la spectra in order to make quantitative statements on 
their potential differences requires the use of such an elaborate 



two-channel deconvolution technique. The code specinholucy 
enables a truly quantitative analysis of high-z SN la spectra 
and opens the path for meaningful detections of potential evo- 
lutionary signatures in their spectra. 

2. The algorithm and its implementation 

2.1. Decomposing the SN from its host galaxy 

iLucv & Walshl ( I2003I) - subsequently LW03 - have described 
two iterative techniques for decomposing a long-slit spectrum 
into spectra of designated point sources and an underlying 
back ground. The methods are based on a two-channel restora- 
tion faook & Lucvlll994t lLucvlll994 using the Richardson- 
Lucy IRichardsoril l 1972: Lucvl l 19741) iterative restoration. The 
essence of this technique is to treat the point source(s) and the 
extended background as two channels for restoration. It must 
be emphasized at the outset that this restoration occurs purely 
in the spatial direction and no implied spectral restoration is 
undertaken. Thus, for instance, residuals from fringing correc- 
tions should not be relevant unless the fringes are tilted with 
respect to the dispersion axis, in which case an additional inho- 
mogeneous component would be added to the background. The 
first channel contains the point source(s) which are restored to 
delta functions using an appropriate PSF. For a spectrum the 
PSF must be specified as a fun ction of wavelength; th is is sim- 
ply the aforementioned SSF in lLucv & Walshl (12003). In order 
to iteratively restore the extended component, it is necessary to 
impose a limiting resolution, larger than that set by the PSF, in 
order to prevent the second channel from modelling the point 
sources as peaks in the extended component. The resolution 
limit for the background is set by defining the restored back- 
ground to be the convolution of a non-negative auxiliary func- 
tion with a wavelength-independent resolution kernel R. Then 
no feature whose width is less than R can appear in the re- 
stored background. If the width of R is greater than the PSF, 
then the convolution of the model of the extended background 
iff with the PSF cannot fit a point source, which must therefore 
be modelled by the first channel. 

Two methods are described in LW03: one for a homoge- 
neous background (i.e. the spectrum of the background is ho- 
mogeneous - does not vary with position); the second, more 
general, case allows the spectrum of the background to vary 
as a function of position along the slit - termed the inhomo- 
geneous background case. In the case of a galaxy including a 
SN the homogeneous case is appropriate. The observed spec- 
trum <j)(A, y), where A is the dispersion direction and y the cross- 
dispersion direction, can then be modelled as: 



-I 



F(A, i])P(A,y- ri)dtl + f SN (A)P(A,y- y w ), 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, operated by the Association of Universities for 
Research in Astronomy, Inc., under contract to the National Science 
Foundation of the United States. 



where 77 is the independent variable for the spatial dimension 
(cross dispersion direction). The first term represents the spec- 
trum of the background and the second represents the spectrum 
of the SN. In this study it is /sn^X the spectrum of the SN, 
which is the desired output product. P(A, y - 77) is the SSF, i.e. a 
PSF as a function of wavelength. The restored background, i//, 
is a convolution of an non-negative auxilliary function x(A, 
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and the non-negative normalized resolution kernel centred at 
i; = f,K(i7-0: 



F(A,t,) = J 

which is convolved with the SSF to match the observed ex- 
tended component. R is specified as a Gaussian function of 
sigma cr kernel . 

The unknowns /sn and y are determined by iterative im- 
provement of the fit of (p to the observed spectrum, using the R- 
L algorithm to improve both the point source spectrum and the 
background until convergence is reached. The iterations steps 
are described in detail in LW03. There are two convergence 
criteria: the fractional change in the spectrum per R-L iteration 
summed for all wavelengths; the fractional change in^ per iter- 
ation summed over all wavelengths and cross-dispersion range. 
In practice, depending on the particular details of the spectrum, 
one of these criteria may converge much faster than the other, 
though both must be met in order for the code to converge. 

2.2. Practical implementation for SN and host galaxy 

The situation of a SN observed in a host galaxy is an ideal use 
of the inhomogeneous case (see LW03 for some other exam- 
ples). This is the simplest application since there is ideally one 
point source and a structurally well-resolved galaxy. When the 
SN is in the outer regions of an early-type galaxy, then the so- 
phistication of the technique is probably not essential and more 
traditional methods, such as a linear or polynomial fit to the 
background and extraction, with or without weighting, of the 
point source is adequate. However, when the SN is near the 
centre of a galaxy, or the galaxy does not have a smooth radial 
profile, or there is line emission in the galaxy host (e.g. for late- 
type galaxies), then the simple methods fail and a dedicated ex- 
traction or a restoration approach is mandated. In section0we 
compare the technique presented here with other less advanced 
methods. 

The inhomogeneous case is implemented in the IRAF code 
specinholucy, which is available in the ST-ECF layered pack- 
age specres. In addition to the input long-slit 2D spectrum, 
there are two important input parameters: the SSF, which obvi- 
ously must match as well as possible the cross-dispersion pro- 
file as a function of wavelength for the SN observation and the 
position of the SN on the slit. The inhomogeneous decompo- 
sition technique requires the position of the point source that 
is to be extracted from the extended background (i.e. the host 
galaxy). Errors in this position will result in mixing of the two 
channels in the output-restored spectrum so that for example 
the background is excavated asymmetrically around the point 
source. However, this requirement is not as limiting as it sounds 
since one of the outputs is a 2D restored version of the input 2D 
long-slit spectrum (including both the SN and the host galaxy) 
which can be directly compared with the latter spectrum and a 
mismatch of the position shows up as a higher frequency com- 
ponent than the resolution kernel in the neighbourhood of the 
designated SN. The primary output is the spectrum of the SN, 
which is strictly ID. If an error image is available for the input 



long-slit spectrum, it can be supplied to produce Monte Carlo 
error estimates for the restored output spectra (point source and 
extended background). 

3. The importance of the slit spread function 

LW03 discuss some of the difficulties with using an empiri- 
cal SSF. For ground-based data a spectrum of a bright standard 
star is ideal, but this must be taken in identical conditions to the 
SN spectrum. In practice this is almost impossible to achieve. 
Even a star on the slit of the same observation as the SN can- 
not be guaranteed to provide an ideal SSF, since the star may 
not be centred on the slit as the SN is centered and instrumen- 
tal optical aberrations can produce off-axis distortions of the 
spatial profile of a point source. A simple task, specpsf, is 
provided to construct model SSFs in the specres package. A 
number of PSFs (viz stellar images) are provided as input at 
specified wavelengths and these are sampled by a slit of spec- 
ified size (relative to the PSFs) and the signal within a slit in- 
tegrated across the slit width as a function of offset position 
(cross-dispersion). The set of point sources are then spline in- 
terpolated in the dispersion direction to provide the 2D spec- 
trum of the point source, the SSF. For an observation with no 
accompanying SSF, a considerable amount of trial-and-error is 
required to choose the seeing of a set of model PSFs in order 
to provide a high quality restoration, and thus successful ex- 
traction of the SN spectrum. In the low signal-to-noise cases, it 
will be the major source of uncertainty on the SN spectrum (see 
section l5Jt . However, all methods which rely on provision of 
a PSF (in imaging) or an SSF (in spectroscopy) are subject to 
this difficulty; it usually provides the fundamental limit on the 
uncertainty of the restoration. 



3.1. Measuring the seeing 

To successfully extract the point source spectrum from the ex- 
tended background component one needs in principle to know 
the width of the point source spatial profile at each dispersion 
coordinate - i.e. the FWHM-wavelength relation for that ob- 
ject. In practice however we often cannot do this due to poor 
S/N (this is especially the case in high-z SN la spectroscopy) 
and so we rely on average quantities, such as the mean seeing 
during the observation. One can then imagine reconstructing a 
more reliable description of the seeing (6) variation along the 
dispersion direction (A) using for insta nce the following seeing- 
wavelength relation dSchroederll987l) : 



= A/r 
r Q oc A 6 ' 5 



<ocA 



-0.2 



where ro is the so-called Fried parameter, which describes the 
quality o f a wave that has p ropagated through atmospheric tur- 
bulence (Fried ll965lfl975l) . 

However, a specific telescope and instrument combination 
can introduce systematic errors due to instrumental spatial dis- 
tortions which are not taken into account in the above rela- 
tion. Fig. [0 shows that one should not rely on measurements 
by external seeing monitors on observatory sites - such as 
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4000 5000 6000 7000 8000 9000 4000 5000 6000 7000 8000 9000 

Observed Wavelength [A] Observed Wavelength [A] 



Fig. 1. Top panel: FWHM of the cross-dispersion profile of 
the standard star LTT 7987 taken with VLT+FORS1 on UT 1 1 
October 2002. The single data point with an error bar shows 
the mean seeing as measured by the DIMM station - namely 
1.55+0.12" (the pixel scale of FORS1 is 0.2" pix" 1 ) at 5000A. 
Overplotted is the FWHM-wavelength relation fixed at the 
DIMM data point, both for the predicted exponent (-0.20) and 
our best-fit value (-0.29). Lower panel: Fractional residuals 
with respect to the 9 oc A~ l relation, for I = 0.20 (open squares) 
and / = 0.29 (filled squares). 

the Differential Image Motion Monitor (DIMM) on the ESO- 
Paranal site, which m easures the image quality via a differential 
centroiding method (Sarazin & RoddieJfl990l ISandrock et alJ 
l2000l) . Although precise, these measurements are inadequate 
for our purpose since the DIMM probes a region of sky signif- 
icantly diffe rent from that where the telescope is pointing (see 
IPatatl 1200.3. Fig. 17). Moreover the optical train is different 
and the observation wavelength is restricted to 5000A. Even 
varying the exponent in the above seeing-wavelength relation 
does not enable to accurately reproduce the FWHM profile of 
a bright standard star (Fig.^. We will see in section l5~3l that 
errors of > 15% in the determination of the FWHM of the point 
source can lead to significant errors in the restored point source 
flux. One possible alternative is to construct a synthetic SSF 
based on the wavelength-dependent FWHM characteristics of 
the point source spectrum to extract. This has the prime advan- 
tage of extracting the point source with a PSF spectrum that 
has been subject to the same seeing variations. To do so one 
can either use the SN spectrum itself or carefully align the slit 
to include a bright single star as well as the object of interest. 

3.2. Generating synthetic PSF spectra 

The IRAF-implemented code specpsf is illustrated here in 
Fig. [2] A standard star is used to determine the correspond- 
ing FWHM-wavelength relation, from which gaussian PSF im- 
ages are generated every 50A. The resulting stellar images are 
run through specpsf and the profile of the corresponding output 



Fig. 2. Illustration of SSF synthesis using specpsf. The input 
spectrum is the same as in Fig.^ Top panel: FWHM of LTT 
7987 (open squares) and of the synthetic SSF (filled triangles). 
The dashed line shows a polynomial fit to the synthetic SSF. 
Lower panel: Residual plots of the FWHM of the synthesized 
SSF (open squares) and polynomial fit to the synthetic SSF 
(filled circles) with respect to the input SSF. 

PSF spectrum (or SSF) is compared with that of the input stan- 
dard star. The errors fall below 0.1% over the whole wavelength 
range. We also show a polynomial fit to the data, since in prac- 
tice the scatter is significantly larger for low S/N point sources, 
and such a fit to the data is more representative of the actual 
seeing variations in the dispersion direction. Here the residuals 
are as high as ~2% in the red region of the spectrum (mainly 
due to an increase of sky brightness in this spectral range), still 
an order of magnitude below the lower threshold where errors 
on the FWHM start to matter (see below). 

3.3. Impact of the SSF on spectral restoration 

The astronomer wishing to extract point source spectra using 
specinholucy has various options for the SSF. In low S/N cases 
like the ones we describe in this paper the choice of the SSF is 
by far the most limiting factor and will have a severe impact on 
the quality of the restoration. 

We illustrate the impact of the SSF choice on the restoration 
of SN 2002go, a Type la supernova at z = 0.236 (IAUC 7994) 
slightly offset from the center of its host galaxy (see also sec- 
tion|6). We consider three options for the SSF: (1) a synthetic 
SSF based on FWHM measurements of SN 2002go itself, (2) 
a PSF star that happened by chance to be on the same slit as 
SN 2002go and (3) a VLT+FORS1 standard star for which the 
FWHM-wavelength relation was the closest we could find to 
that of SN 2002go. The resulting cross-dispersion profiles are 
shown in Fig. [3] and the corresponding mean seeing values are 
shown in Table [2 Note that the quoted mean values are not 
representative of the wavelength-dependent nature of the see- 
ing and one should produce plots such as those shown in Fig.|5] 
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Fig. 3. Examples of SSFs used to extract SN 2002go: (1) the 
supernova itself; (2) a bright PSF star that happened to be on 
the same slit; (3) the closest FORS 1 standard star to match the 
supernova profile. Top panel: The dashed line is a polynomial 
fit to the SN profile, from which low S/N points affected by the 
increase in sky background have been excluded. Lower panel: 
Fractional differences in the SSF FWHM with respect to the 
polynomial fit. 

Table 1. SSF FWHM data 



SSF 


pixels" 




arcseconds" 


SN 2002go 


o 014.74 


+ 0.41 


77 M5 

' '0.70 


±0.08 


star on slit 


4 7 1 544 


±0.39 


94110 


±0.08 


standard star 


A 174.64 

'3.93 


+ 0.23 




±0.05 



(a) the quoted values are read as < mean > max ±lcr error 

v ' 1 mill 



to appreciate the differences between the SSFs. Despite having 
been observed on the same location of the CCD chip, the SN 
and the standard star were not observed on the same night and 
were thus subject to different seeing conditions. The difference 
is also a function of wavelength, hinting at possible instrumen- 
tal effects on the spatial profile. These same comments apply 
to standard stars having been observed on the same night as the 
object of interest, and we discourage the potential user of this 
method from using standard stars for the SSF. 

More interesting, perhaps, are the quasi-systematic differ- 
ences in the profiles of the SN and the PSF star that was acci- 
dentally placed on the same slit. The SN and PSF star should 
be subject to identical seeing variations, but the reconstructed 
profiles differ by as much as ~ 30%, or 0.2". There could be 
several reasons for this: (a) the observation is not made at the 
parallactic angle, and so the dispersion direction is not along 
the slit; (b) the star spectrum is located at a different location 
on the CCD chip where the spatial distortion differs; (c) the 
star is not centred on the slit and we are not measuring its 
true PSF profile. FORS1 is equipped with a longitudinal at- 
mospheric dispersion compensator (LADC) which has an ac- 
curacy < 0.1" (0.5 pix) over the considered wavelength range 



jAvila etalJl997t) . and so (a) should not be relevant. FORS 1 is 
linear at the < 0.4% level (see the FORS 1+2 User Manual 2 ) 
and distortions associated with the instrument's optics should 
not affect the whole spectral range at this level, so (b) should 
not matter either. This leaves (c) as the probable explanation 
for this difference, which is significant enough to be a major 
source of error in the restoration. 

The impact of these various SSFs on the restored SN flux 
is shown in Fig. [4] We view the solid curve as the most accu- 
rate restoration, and one sees that using the other two (wider) 
SSFs means we are also restoring part of the extended back- 
ground component (sky and host galaxy) in the point source 
spectrum. In the case of the star that was by chance located 
on the same slit the ~ 20% difference in the FWHM results 
in an error in the restored SN flux that can be as high as 100- 
150%! We see from the two inlays in Fig.|4]that the slope of the 
pseudo-continuum around strong SN la spectral features (Ca II 
H&K and Si II) changes significantly due to contamination of 
the underlying background spectrum, and so the very defini- 
tion of "line strength" in this case is very ambiguous. The fact 
that the FWHM of the point source varies with wavelength will 
cause variable contributions of the background to the restored 
spectrum, and will affect the determination of line ratios of dis- 
tantly spaced lines, if the PSF vector is not determined prop- 
erly. Thus, the extraction has a significant impact on empirical 
correlations based on such measurements, as the errors made 
at this stage will add to those associated with the subsequent 
calibratio n in flux. Examples of such correlations are given in 
iNueent et all 19951 and are used to determine distances to SNe 
la in lRiess etalJl998bl 

One might ask how one can decide which SSF to choose 
to restore a particular point source spectrum and if it is pos- 
sible to a posteriori tell whether the chosen SSF was indeed 
the appropriate one. Indeed, through restoring the complete 2D 
background spectrum specinholucy enables us to check the ac- 
curacy of the restoration at each pixel of the input 2D frame, 
or simply in the spatial (cross-dispersion) direction by compar- 
ing the collapsed spatial profiles of the input and restored 2D 
spectra. This is illustrated using real data in Fig. 1121 

4. The spatial resolution kernel 

The choice of the width of the resolution kernel for the back- 
ground resides with the user. It is implemented as a Gaussian 
of user-chosen sigma cr kernel in the specinholucy code. It is dif- 
ficult to give a general guide to its choice. Obviously it should 
be wider than the equivalent maximum Gaussian sigma of the 
SSF crssF.max in the cross-dispersion direction. The sigma of 
the kernel has to be tuned and is very dependent on the nature 
of the extended source and more specifically on the extended 
background component in the vicinity of the point source. 
Comparison of the restored background with the input spec- 
trum will quickly show if narrow features in the extended back- 
ground, above the noise, have failed to be modelled through 
choice of too wide a kernel width (see Fig. [5] left column). If 
the kernel width is too narrow then the point source will be 

2 http://www.eso.org/instruments/fors 1/ 
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Fig. 4. Impact of the SSF on the restored point source spectrum 
of SN 2002go. Refer to Fig.[3]for the meaning of the line anno- 
tations (1), (2) and (3), noting that (1) corresponds to the poly- 
nomial fit to the SN 2002go profile. The two inlays zoom in on 
the blueshifted Ca II H&K and Si II features, both prominent in 
SN la spectra, and highlight the impact of non-optimal spectral 
extractions on empirical correlations involving line strengths 
(see text). 
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Fig. 5. Impact of the width of the spatial resolution kernel 
enamel on the background restoration. Top panel: Input aver- 
age (dashed line) and restored background (solid line) cross- 
dispersion profiles. Middle panel: Residuals of the previous 
plot, revealing the point source channel (i.e. the supernova it- 
self). Lower panel: Spatial residuals of the input and restored 
2D frames. Negative (positive) values in the residual plots 
mean we are over (under)-restoring the flux at that location. 



partially modelled as a peak in the background channel and 
the fluxes in our restored point source spectrum will be un- 
derestimated (Fig. |3 right column). The optimal width for the 
spatial resolution kernel corresponds to a compromise between 
the lower spatial frequency of the extended background cross- 
dispersion profile and the higher spatial frequency of the point 
source cross-dispersion profile. A decent first guess would be 
Ckemei ~ 2-3<xssf- I n Fig- ID (middle column) we have used 
0-kemei ~ 2.3cr SSF > <T S sF,max- From Fig.|3]we further note that, 
if one can clearly see the impact of a wide kernel on the spatial 
residuals (5 ga i+SN, discriminating between the narrow and opti- 
mal kernel widths - respectively 1.5 and 3.8 - is not obvious 
a priori. The optimal kernel width in Fig. [5] was chosen such 
as to maximise the counts in the point source channel, or 5 ga i, 
whilst keeping the spatial residuals 5 ga i+SN below the statistical 
noise limit of the input 2D frame. 

One important test is to ensure we are not capable of restor- 
ing a PSF-like spectrum from a pure extended source. This 
is particularly relevant to high-z SN la spectroscopy since in 
many cases one is only able to extract a spectrum of the su- 
pernova host galaxy whereas one is convinced of having ob- 
tained a combined spectrum including the SN. Then the cur- 
rently widespread technique of subtracting a galaxy template 
from this "combined" (and noisy) spectrum can sometimes re- 
veal SN-like features where no supernova is present! This is 
illustrated later on in section l73l where we compare the specin- 
holucy output with that of a statistical algorithm {SN -fit) which 
decomposes a ID spectrum into its SN and galaxy components 
using a galaxy spectral template (see also Fig. 1141 . 



Here we have used the uncontaminated galaxy spectrum of 
the host of SN 2002bo (see section obtained by construct- 
ing a mirror image of the portion of the galaxy profile devoid 
of supernova signal. We try to extract a point source compo- 
nent from this pure galaxy spectrum at different locations of 
the galaxy trace, including the centroid, using for this purpose 
three different SSFs of widths corresponding to the mean and 
extrema seeing values measured by the DIMM station for this 
observation. We then compare the restored point source flux 
with the statistical noise limit of the input 2D spectrum at the 
position of the point source. Should this ratio fall below one 
this means that no point source was detected. We see from Fig. 
|6]that in all cases the restored spectrum is inconsistent with that 
of a PSF-like object over the whole spectral range. 

5. Testing the algorithm on simulated data 

In this section we present tests of the algorithm on simulated 
data, namely 2D spectra of a Type la supernova (SN la) at 
Z = 0.5 embedded in a late-type galaxy. We have chosen 
to reproduce a situation where a one hour-long exposure is 
taken at the ESO Very Large Telescope VLT-UT1 (8.2m) us- 
ing the FOcal Reducer Spectrograph (FORS1) in Long-Slit 
Spectroscopy (LSS) mode. The grism corresponds to 300V, and 
the slit width is 1". The reason behind chosing these specific 
settings is that they correspond to those used for the real data on 
which the algorithm will be tested further in section |6] In what 
follows we assume a (Q.m, Qa, h) = (0.3, 0.7, 0.72) cosmology 
(where h = ■—- — ^t— — r is the dimensionless Hubble con- 

v 100 km s 'Mpc 1 

stant), and so the luminosity distance corresponding to z = 0.5 
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Fig. 6. Restored point source spectra in units of the statistical 
noise of the input pure galaxy spectrum. The different plots 
correspond to different positions of the point source to extract, 
marked by a filled square in the galaxy spatial profiles shown 
in the right column. The different lines correspond to different 
seeing conditions. The horizontal dashed line marks the limit 
below which no point source was detected in the background. 
Note the change of scale in the lower plot, where no flux at all 
is restored in many bandpasses. The peak above the dashed line 
in the uppermost plot is due to bad cosmic ray removal. 

Table 2. Distances and flux scaling 



galaxy d (Mpc) Reference factor 2 
NGC 4526 13 Drenkhahn & Richtler 1999 1.5 x lfr 5 
NGC6181 34 Sil'chenko et al. 1997 1.1 x 10~ 4 

(a) Factor used to scale down the flux, computed as 

W4, : =0.5) 2 X(l+ Z r'. 



5.1 .1 . Supernova spectrum 

The 2D supernova spectra are synthesized using the ID, flux- 
calibrated spectra of the Type la supernova SN 1994D in NGC 
4526 (cz = 448 km s -1 ) as obtained from the SUSPECT 
database 3 . Where spectral coverage at a given phase was lack- 
ing (i.e. at -6d, +0d, +6d, +8d and +14d), approximate syn- 
thetic spectra were obtained from the two observed spectra 
closest in phase using the UBVRI photometry published in 
Pata t et al.ldl996l) . The same technique was applied to observed 
spectra whose wavelength range did not extend below 3650A 
in the blue (at -8d and +10d), so as to make sure to cover the 
Ca II H&K (3934A,3968A) absorption trough bl ueshifted to 
~375 0A, prominent in SNe la optical spectra (see iFilippenkol 
1 19971 for a review). Each spectrum is then redshifted to z — 0.5 
and rebinned to the pixel scale of FORS 1 equipped with a 300V 
grism and a 1 "slit (~ 2.66A pix -1 ). The flux (in units of erg s" 1 
cm -2 A -1 ) is scaled down to account for the change in luminos- 
ity distance as the source is artificially placed at z = 0.5 - as- 
suming SN 1994D to be at the same distance as its host galaxy 
- and is converted to counts (ADU) using a FORS1 sensitiv- 
ity function. At this stage we have a ID supernova spectrum 
in ADUs as a function of wavelength, corresponding to a one- 
hour integration on FORS 1 . Finally, a 2D supernova spectrum 
is synthesized by replicating this ID spectrum along the lines 
of a 2D image array (i.e. along the dispersion axis) and multi- 
plying each image column by a ID gaussian profile normalised 
to unity of FWHM corresponding to 1 "seeing. Note that the 
flux vector is normalised to an airmass of 1, and no scaling was 
applied to account for observations (usually) made at higher 
airmasses. 

The supernova data used in this simulation are summarised 
in Table [3] The quoted /?-band magnitudes are only approxi- 
mate and are affected by systematic effects due to the lack of 
any ^-correction. For comparison the mean number of counts 
(in ADUs) per pixel is 105.7 and 19.8 for galaxy and sky (in- 
cluding sky lines), respectively. 



is c/l,z=o.5 ~ 2.72Gpc. The factors used to scale down the fluxes 
as the SN and galaxy are artificially redshifted to z — 0.5 are 
listed in Table 13 The additional factor of (1 + z) comes from 
the fact that we are integrating the flux over a finite bandpass, 
where the wavelength axis has been diluted by that same factor. 

5.1. Simulated data 

The simulated 2D spectra are a combination of a 2D supernova 
spectrum and a 2D background spectrum, itself consisting of a 
galaxy and a sky spectrum. We vary the phase of the SN, i.e. 
the brightness and the input spectrum, and its position within 
the galaxy from 0.75 to 1.75, in units of FWHM of the galaxy 
trace. Random poisson noise is added to the image using the 
gain and read noise characteristics of FORS 1 so as to degrade 
the overall S/N of the image and thereby reproduce plausible 
observing conditions. 



5.1.2. Background spectrum 

For the galaxy spectrum we choose to artificially place SN 
1994D in NGC 6181, an SC ga laxy part of Ken nicutt's spec- 
trophotometric atlas of galaxies dKennicuti 1 992l) . The reason 
behind this is the prominent [O II] emission at 3727 A present in 
its spectrum (see Fig. [7}, which could then serve as an diagnos- 
tic tool for probing the efficiency of our algorithm in restoring 
a clean point source spectrum. If successful, the restored ID 
spectrum of the point source should be devoid of [O II] flux 
residuals. 

The spectrum of NGC 6181 was redshifted to z = 0.5 and 
rebinned to the pixel scale of FORS 1 using the same procedure 
as for the supernova spectrum. The flux was scaled down to 
account for the change in luminosity distance, and converted to 
counts (again in ADUs) using the same sensitivity function as 
for the supernova. The resulting galaxy has a magnitude m g - d \ ~ 
21.5 in a reshifted B-band filter (very roughly corresponding 

3 http://bruford.nhn.ou.edu/~suspect/. 
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Table 3. SN 1994D data used in the simulation 



phase a 


date b 


m B c 


"^,c=0.5 d 


counts e 


-10 


11/03/94 


13.26 


24.86 


13.5 


-8 


13/03/94 


12.67 


24.27 


22.7 


-6 




12.26 


23.86 


30.7 


-4 


17/03/94 


11.99 


23.59 


37.4 


-2 


19/03/94 


11.86 


23.46 


39.5 


+0 




11.84' 


23.44 


51.6 


+2 


23/03/94 


11.89 


23.49 


48.4 


+4 


25/03/94 


11.99 


23.59 


34.8 


+6 




12.13 


23.73 


31.3 


+ 8 




12.33 


23.93 


27.9 


+ 10 


31/03/94 


12.55 1 " 


24.15 


21.5 


+ 12 


02/04/94 


12.80 


24.40 


15.0 


+14 




13.10 


24.70 


12.7 



(a) SN phase in rest-frame days from B-band maximum. 

(b) UT date of spectroscopic observation. Absence of date means the 
spectrum is synthetic. 

(c) S-band magnitude taken from Table 1 in lPatat et alJll993) . Data 
marked with T results from a polynomial fit to the 6-band light curve. 

(d) the S-band at z = 0.5 is close to rest-frame S-band, and so 

m R (z = 0.5) ~ m B - 2.5 log (^\ * m B + 1 1.60. 

(e) mean total counts in the whole SN profile in ADUs per pixel 
along the dispersion axis. 




distan ce r from the nucleus Jde Vaucouleursl[l959t iFreemanl 
One can derive a global analytical form of the surface 
brightness profile (r) of spiral galaxies, as a function of the 
effective radius r e of the bulge only. In doi ng so we assume 
abulge- to- total luminosity ratio B/T = 0.4 jRatnatunga et alJ 
1999) and a ratio between the ef fective radi i of the bulge and 
disk c omponents r^ e /r^ e = 0.5 JKentHl985llRatnatunga et all 
1999): 



S 5 (r) = Zfc(r) + Zrf(r) 

= 0.76931 S s , e exp 

+2.9343 I. s ,e exp 



-7.6692 
1.3945r 



1.6617r 



1/4 



5500 6000 6500 7000 7500 8000 
Observed Wavelength [A] 



8500 



where Xf,(r) and 2</(r) are the surface brightness profiles of the 
bulge and disk components, respectively, and S s e (r) is the sur- 
face brightness at the effective radius r e . 

The angular size of NGC 6181 is 2.5'xl.l' (NASA/IPAC 
Extragalactic Database 4 ), which corresponds to 4.5"x2.2" at 
z = 0.5 , and in turn to 22.5x11 pixels on a FORS1 image. 
Assuming the slit was placed along the major axis of NGC 
6181 we can assume that the effective radius of the disk r e is 
22.5 pixels. The derived luminosity profile is then normalised 
to unity and is multiplied into each column of the 2D galaxy 
image array. Each column is further convolved with a gaus- 
sian profile of FWHM corresponding to 1" seeing to reproduce 
similar observing conditions as for the supernova. The result- 
ing spatial profile / ga i(r) is shown in Fig. [8] Notice how at these 
redshifts the core spatial profile of galaxies is entirely seeing- 
dominated, and one must rely on the broad extended wings of 
the profile to differentiate between a galaxy and a point source. 

The sky 2D spectrum (the second background component 
for ground-based observations) is generated by uniformly repli- 
cating an observed ID sky spectrum with the required settings 
(300V grism, 1 " slit) along lines of an image array, and scaled 
to the required exposure time (3600s in our case). It is thus by 
construction perfectly flat and monochromatic in the spatial di- 
rection. We do not have to worry about spatial distortions of sky 
lines and sky background removal in the input 2D spectrum is 
rendered trivial. 



Fig. 7. Absolute flux-calibrated galaxy and sky spectra used in 
the simulation. The spectra are shown in observed wavelength 
at z = 0.5, and without shot noise for sake of clarity. The promi- 
nent [O II] emission and the Ca II H&K absorption in NGC 
6181 are labelled accordingly. Overplotted (dashed lines) are 
the spectra of SN 1994D at maximum intensity (top) and ten 
days before maximum (bottom). 

to the S-band at z — 0.5) i.e. more than an order of magnitude 
brighter than when the supernova is at its faintest (see Tables- 
Next, a galaxy spatial luminosity profile was constructed 
based on a standard bulge+disk model. The bulge component 
is usually described analytically using the so-called r 1 ^ 4 law 
and the disk component follows an exponential decline with 



5.2. Simulation steps 

We now move on to outlining the steps made in carrying out 
the simulation. Far from being redundant these are the same 
steps that one should in principle go through when applying 
this technique to real data. They include determining the posi- 
tion (and slope) of the supernova spectral trace, its FWHM as a 
function of wavelength, synthesizing a corresponding SSF and 
running the specinholucy code with the optimal settings. For 
the purpose of this simulation and for ease in computing statis- 
tics these steps have been made automatic, though we strongly 
advise to make them highly interactive when dealing with real 
data. 

4 http://nedwww.ipac.caltech.edu/ 



S. Blondin et al.: Extracting clean supernova spectra 



9 



Profile wings 




0.5 1.0 1.5 

Distance from galaxy centroid [FWHM 



2.0 



gal 



Fig. 8. Smoothed galaxy spatial profile used in the simulation, 
resulting from the convolution of a standard spiral galaxy lu- 
minosity profile and a 1" seeing gaussian kernel. The abscissa 
is in units of FWHM of the convolved galaxy profile (= 3.06 
pix). The effective radius r e in these units is at ~ 7.36. Filled 
squares indicate the position of the SN used in this simulation, 
namely 0.75, 1.00, 1.25, 1.50 and 1.75. 



5.2.1 . SN trace position and slope 

One of the inputs to the specinholucy routine is a position table 
containing the (x, y,slope pix -1 ) of the supernova spectral trace. 
Since in the case of high-z SNe la we are limited in S/N we 
cannot simply give an approximate position and rely on cross- 
correlations with the SSF to determine a more accurate one. 
Here we must input the exact SN trace position to which the 
SSF profiles will be shifted. The results of our simulation (see 
section l5~3l show the impact of positional accuracy on the point 
source restoration. 

To some extent one should base the method used to de- 
termine the trace position on the characteristics of the instru- 
ment used for the observation. In the case of FORS 1 we know 
from observations of standard stars that the position of a point 
source does not vary by more than one pixel across the whole 
CCD chip. Thus, a ~ lpix accuracy can almost routinely be 
achieved by summing up each of the spatial channels and find- 
ing the maximum intensity (other than the galaxy) in a given 
cross-dispersion region. Clearly this will be problematic for 
cases where the SN is close to the centroid of the galaxy 
trace or when its phase is far from maximum (-10d or +14d). 
Alternatively one can try to fit the background around the SN 
trace and subtract this fit from the overall profile. The residuals 
can then be fit with a gaussian profile and the trace position is 
taken to be the profile centre. 

The result from this operation is a relation between SN 
trace position and dispersion coordinate x, which is fit linearly 
so as to determine the trace position and slope. We now have 
our input (jc, }',slope pix -1 ) table necessary to run specinholucy. 
A plot of the error made in determining the trace position using 
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Fig. 9. SN trace position residuals as a function of the SN 
phase, for different locations of the SN with respect to the 
galaxy centroid. Each curve corresponds to a specific SN posi- 
tion, given in units of FWHM ga |. 

this approach is shown in Fig. [9] The position is best deter- 
mined when the signal-to-noise ratio is high, namely when the 
contrast is high between the SN and the underlying background 
, as expected. 

5.2.2. SN trace FWHM and SSF synthesis 

To construct a synthetic SSF we need to know the FWHM of 
the SN trace as a function of the dipersion coordinate. We pro- 
ceed in a way analogous to determining the SN trace position, 
except we are now interested in the width of the gaussian fit 
to the residual signal left over after sky and galaxy subtrac- 
tion rather than the gaussian centre. A similar relation between 
FWHMsn and dispersion coordinate is established and a table 
of (FWHMsn, wavelength) values is elaborated and serves as 
input for the SSF synthesis using specpsf. A plot of the error 
made in determining the trace FWHM using this approach is 
shown in Fig.^| Again, the FWHM is better determined for 
cases with a higher SN flux relative to the galaxy background. 

5.2.3. Running specinholucy 

Once we have determined the position (and slope) of the SN 
trace and synthesized an SSF spectrum using its FWHM- 
wavelength dependency, we can run specinholucy on the 
inpu t sky-subt r acted 2D spectrum. We refer the reader to 
iLucv & Waishl d2003h and to the IRAF help pages for the 
specinholucy task for a more complete description of the input 
parameters. A centroid fitting is used to determine the SSF peak 
at each column, and no subsampling (in the spatial direction) 
is performed since the resulting PSF profile is not oversampled 
with respect to the input 2D spectrum. A cubic spline is used 
as the interpolation method to shift the SSF to the position of 
the point source - which is defined exactly in the input posi- 
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Fig. 10. SN trace FWHM residuals as a function of the SN 
phase, for different locations of the SN with respect to the 
galaxy centroid. The curves and plotting symbols have the 
same meaning as in Fig. [9] 



tion table. Note that no input statistical image to determine the 
statistical error at each pixel is used in these simulations. 

The only input parameter made to vary from one run to 
the other is the width <x kerne i of the gaussian used for the spa- 
tial resolution kernel. We see from Fig.[8]that the galaxy spatial 
profile is entirely seeing-dominated, meaning that the widths of 
the SN and galaxy traces are comparable at each wavelength. 
This in turn implies that the spatial frequency of the galaxy pro- 
file will vary hugely when varying the SN position from 0.75 
to 1.75 FWHMg a i. To avoid confusion in the point/extended 
source discrimination we must use a narrow (wide) spatial res- 
olution kernel when the SN is close to (far from) the galaxy 
centroid, where the spatial frequency is high (low). The input 
values for the smoothing kernel were chosen such as to min- 
imise the mean flux residuals in restoring all phases of the SN 
spectrum at a given position. For all runs the convergence cri- 
teria for the SN and background spectra per R-L iteration (see 
section|2j were set to 0.1% and 1%, respectively. 

Two runs (A & B) were executed: in run A we use the SN 
position and FWHM as determined using the method outlined 
in sections r5.2.1l and l5.2.2l whilst in run B we use the known 
values of the SN trace position and FWHM to synthesize the 
SSF and restore the SN spectrum. This enables us to investigate 
the joint impact of the errors on the SN position and FWHM on 
the accuracy of the SN spectral restoration (see section l5~3l . 

5.2.4. Statistical calculation 

To evaluate the efficiency of the method we compare the 
(known) input SN spectrum with its restored version. More 
specifically we compute the ratio of the flux residuals with the 
statistical variation in the input 2D spectrum. Should this ratio 
fall below one this means that we are statistical noise-limited 
and cannot improve on the restoration. 



In this simulation we compare the flux residuals both over 
the whole observed spectral range (5500A-8500A) and around 
the region of host galaxy [O II] 3727 A emission. This region 
is simply defined in wavelength space (in A) as 3727(1 + z) + 
FWHM [0 ii] (1 + z), where FWHM [0 n] is the width of the [O 
II] line in the input non-redshifted spectrum of NGC 6181 and 
z = 0.5 is the redshift of the simulation. For both cases the flux 
residual 6F is evaluated at each pixel according to: 



5F 



|f SN,i - Fst. 



V^SN.i + -Fgal.i + -Fsky.i 



where FsN.i and are the input and restored SN fluxes, re- 
spectively. F gai,; and F s k y ,i are the galaxy and sky fluxes at the 
location of the SN, defined as an interval centred on the exact 
SN trace position and of width the FWHM of the SN trace. It 
is important to compute the statistical noise as that in the to- 
tal signal at the position of the supernova (i.e. including the 
underlying background) and not just in the SN signal itself. If 
6F < 1 one can in principle consider the SN restoration to have 
converged. 

5.3. Simulation results and discussion 

The simulation results are summarised in Fig. ^2 Both runs 
A and B (see section l5.2.3l > are shown in order to evaluate the 
impact of the error made on the determination of the SN trace 
position and FWHM on the restoration of the supernova spec- 
trum. 

These plots clearly show that the limiting factor of this two- 
channel restoration technique (as for any other spectral extrac- 
tion method) is the contrast between the point source and the 
underlying background. The overall flux is systematically re- 
stored to the statistical noise limit when the supernova is at its 
brightest stages (-5d < phase < +5d), irrespective of its po- 
sition with respect to the centroid of the galaxy trace. In this 
phase regime the error on the SN trace position and FWHM is 
also at its lowest (see Figs.l9land ll0t and has little impact on 
the quality of the restoration - as shown by the dashed lines in 
Fig. ^2 This is not always the case for the [O II] flux residuals, 
where the added background noise due to the presence of this 
emission line affects the restoration (cf. at positions 0.75 and 
1.00, in units of FWHM gal ). 

One sees the dramatic impact of the error made on the SN 
trace position and FWHM at low S/N, namely when the SN is 
outside the -5d < phase < +5d range and close to the centroid 
of the galaxy (at positions of 0.75, 1.00 and 1.25). For errors 
in the position > 0.3 pix and FWHMsn <; 15%, the restora- 
tion fails to reach the statistical noise limit. These are the main 
sources of systematic error in the restoration method and great 
care should be taken when determining the SN trace position 
and FWHM when applying it to real data. 

6. Testing the algorithm on real data 

In this section we apply this restoration technique to real su- 
pernova data. An outline of the steps involved (determination 
of the SN trace position and FWHM; SSF synthesis; running of 
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Fig. 11. Top panel: Flux residuals 6F as a function of the SN phase for different positions of the SN with respect to the centroid 
of the galaxy trace, given in units of FWHM ga |. 6F < 1 means the SN spectrum has been restored to the statistical noise limit. 
The solid and dashed lines correspond to runs A and B, respectively (see section l5.2.3> . Lower panel: [O II] flux residuals as a 
function of phase for different SN positions. 



specinholucy) has already been given in section l5~2l We simply 
provide the input parameters that were used to run the code 
(see Table |4}. All the data were collected by the authors via 
the ESO observing programme 170.A-0519 and are available 
through the ESO Science Archive Facility 5 . The observations 
were made with VLT+FORS 1 in Long Slit Spectroscopy (LSS) 
mode with a 300V grism and a 1" slit. The results of the two- 
channel restoration are shown in Figs. [21 and ^] whilst the 
individual objects are presented in Table|5] On average ~ 150 
iterations were needed to reach convergence for the restoration 
of the point source spectra, corresponding to ~ 2 minutes on a 
2.4 GHz Pentium 4 processor. 

There is one important distinction to be made in the calcu- 
lation of the flux residuals here as opposed to that made for the 
simulation in section l5".2.4l For observed cases we do not know 
how the total flux F tot separates into the individual components 



5 http://archive.eso.org/ 



F sn, ^ g ai and F s ^ y , and so our evaluation of the accuracy of the 
restoration is restricted to the calculation of the mean total flux 
residuals 6F tot : 

Ftot.i - F tot ,rl 
Or tot = —= . 

V-Ftot.i 

This is not as limiting as it may seem since we are again able 
to quantitatively evaluate our combination of <T) esm& \ and ctssf 
through maximising the counts in the point source channel 
whilst ensuring that 6F tot < 1 (cf. section^. 

Fig.ll2lillustrates what has been said about an a posteriori 
check of the adequacy of both the SSF and the spatial resolu- 
tion kernel. In all cases we not only fit the SN trace but also the 
underlying background such that the flux residuals 6F tot < 1 . 
For the case of SN 2002bo a secondary point source (indicated 
by an arrow) in its vicinity had to be included in the input po- 
sition table - and hence allocated to the point source channel - 
for the restoration to succeed. 
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Table 4. specinholucy parameters 



SN a 


^"kernel 


FWHM 


SN+bg c 


SN/bg d 




[pixels] 


DIMM 


SSF 






2002bo_l 


4.0 


6.25 


6.93 


63.6 


4.8 


2002bo_2 


4.0 


7.95 


6.57 


62.5 


4.6 


2002gr 


2.2 


5.85 


3.46 


9.0 


1.4 


2002go 


3.5 


5.75 


3.97 


5.2 


2.1 



(a) IAU designation. Note that two spectra were taken of SN 2002bo. 

(b) mean FWHM in pixels of the gaussian seeing corresponding to 
the DIMM station measurement or to the SSF used to extract the 
spectrum. 

(c) mean number of counts (ADU) in the dispersion direction of the 
input 2D spectrum at the location of the SN. 

(d) ratio of the integrated SN flux and that of the underlying 
background. One can calculate the mean signal-to-noise per pixel 
using SNR pix" 1 = SN/bg + 1. 



Table 5. Supernova observation summary 



SN a 


Archive id. b 


UT date 


f 


phase d 


2002bo 


2002bo 


06 Dec 2002 


0.0042 


+258 


2002gr 


sloan9 


1 1 Oct 2002 


0.091 


+ 14 1 ' 


2002go 


sloan5 


1 1 Oct 2002 


0.236 


+2+ 



(a) IAU designation. 

(b) target name for search in the ESO archive. 

(c) redshift determined from nebular lines in the host galaxy. 

(d) SN phase in rest-frame days from fi-band maximum. Phases 
marked with + have been determined from cross-correlations with 
local SN la templates using the SuperNova IDentification (SNID) 
code fronrv et all2005l) . 

Fig. [D] on the other hand shows the flux residuals in 
the dispersion direction, and one can immediately appreciate 
the successful allocation of (extended) atmospheric absorption 
features and sky emission lines to the background channel. 
Moreover one can immediately pick out the spectral regions 
affected by systematic errors in the restoration (6F tot > 1) and 
only use those where 6F tot < 1 for analysis. The (spectral) 
residuals in this case were evaluated in the dispersion direc- 
tion by evaluating the spatial residuals in a region centred on 
the SN position and of extent 3cr SSF . This ensures that we are 
only including the point source and its immediately underlying 
background. 

7. Comparison with other methods 

We choose to compare our method with alternative techniques, 
namely a standard extraction using IRAF, an iterative gaussian 
extractor and a statistical decomposition of a ID flux-calibrated 
spectrum into SN and galaxy contributions. The reason for in- 
cluding the latter (ID) approach is to highlight the need for the 
SN-galaxy separation to occur at the earliest stages of the re- 
duction process, and not as part of a post-processing chain of 
data analysis. The results are summarised in Fig.^J Note that 
algorithms similar to ours do exist (though unlike ours are not 
publicly available), and should thes e be of interest to the reader 
we refer him/her to lLucv & Walsh ( 2003), section 6. 



7.1. Standard extraction in IRAF 

For most purposes IRAF and similar reduction software are 
sufficient for extracting a high S/N supernova spectral trace. 
The extraction can further be optimised by e.g. tracing the 
point source signal along the dispersion direction or variance- 
weighting the output flux. However, the hold on systematic er- 
rors is null and one tends to evaluate the resulting spectrum on 
purely qualitative grounds. In the case of high-z SN la spec- 
troscopy the extraction of the SN spectrum is often considered 
successful when the cross-correlation with a local SN la tem- 
plate spectrum is maximal. This clearly bypasses the specificity 
of the high-z spectrum and affects the search for systematic 
spectral differences between high-z SN la spectra and their lo- 
cal counterparts. 

Of equal concern is the arbitrary definition of the SN "back- 
ground" in IRAF. One specifies a set of spatial coordinates in 
the vicinity of the SN trace to be fit by some user-defined poly- 
nomial . . . and essentially hopes for the best! In particular one 
has to find a compromise between the varying width of the SN 
cross-dispersion profile (i.e. its FWHM-wavelength relation) 
and the varying structure of the underlying background. Cases 
such as that encountered for SN 2002bo - namely the con- 
tamination of the close SN background by a secondary point 
source - are extremely delicate to handle in a coherent way 
with standard spectral reduction packages and lead in this case 
to > 40% errors in the SN flux. Other errors in the extraction 
using IRAF are due to the varying spectrum of the underly- 
ing background, to strong sky emission lines and to the overall 
increase in sky noise at long wavelengths. The extracted spec- 
trum of SN 2002go for instance is clearly affected by sky resid- 
uals, and more specifically the Si II feature (6355A) blueshifted 
to ~ 6150A and characteristic of Type la supernovae lies in the 
same region as the atmospheric A-band at z = 0.236, and is 
clearly contaminated by it. Any measure of line strengths and 
ratios in this region will thus also be affected, leading to signif- 
icant errors in empirical correlations based on such mea sure- 
ments, e.g. the "Nugent relation" (see lNugent et alJll995j and 
section l3~3*l in this paper). In cases where the signal-to-noise 
ratio is low (cf. SN 2002gr) there is obvious contamination 
from the underlying background in the SN spectrum, leading 
to significant errors; in fact, no SN flux at all is output in many 
bandpasses, where specinholucy is able to restore the SN con- 
tribution over the whole wavelength range. 

7.2. Gaussextract - an iterative gaussian extractor 

In applying specinholucy to the restoration of supernova spec- 
tra we need to degrade the resolution of the background chan- 
nel to discriminate between the SN and its underlying back- 
ground. To illustrate the advantages of using a two-channel 
restoration we have written a software (gaussextract) that iter- 
atively extracts the point source spectrum using the same syn- 
thetic SSF as for the specinholucy run. As in specinholucy each 
dispersion channel is treated independently and the process is 
iterated till some minimal fractional change in the restored SN 
flux is achieved. For the sake of comparison with specinholucy, 
gaussextract runs the same convergence criteria for the point 
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Fig. 12. Top panel: Normalised average spatial profiles of the input (dashed line) and restored (solid line) 2D background spectra. 
Two runs of specinholucy were executed, one including the point source in the restored 2D background spectrum and the other 
excluding it, so as to appreciate how well the background underneath the SN was fit. Lower panel: Wavelength-averaged spatial 
residuals in units of the statistical noise of the input 2D spectrum. For all cases we have 6F tot < 1, meaning the combined 
SN+background spatial profile is restored to the statistical noise limit over the whole spatial range. 
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Fig. 13. Top panel: Normalised restored point source spectra (solid line) and underlying background, both including (dashed 
line) and excluding (dotted line) the point source. The spectra have been normalised to the integral of the underlying background 
flux. The symbol ffi denotes the atmospheric A-band (~7600-7630A), whilst arrows indicate strong sky emission lines. Lower 
panel: Spectral residuals in units of the statistical noise of the input 2D spectrum. 



source flux. The difference here is that the SN and underlying 
background are jointly extracted in the same channel, leading 
to significant contamination by sky lines and host galaxy in the 
output SN spectrum (cf. SN 2002bo in FigJH}. 

7.3. Statistical approach using SN-fit 

SN-fit is a software developed by Gregory Sainton (see 
ISaintonl Eo04) to rapidly identify the supernova type and red- 



shift in the SuperNova Legacy Survey (SNLS). The goal is to 
quickly identify candidates from their spectrum, for which the 
SN phase and the fraction of host galaxy light can be estimated. 
The algorithm is based on a x 1 fitting of the observed (and re- 
duced) ID spectrum to known template spectra. More specif- 
ically the observed spectrum is compared with a model made 
of a fraction a of SN (S) and a fraction /3 of galaxy (@). The 
following model spectra M are built : 
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Fig. 14. Top panel: Restored supernova spectra using specinholucy (solid line) and other methods presented in sectional The 
spectra are normalised to the maximum flux value of the specinholucy output. Due to the specificities of SN-fit - namely the 
restriction of the local SN template spectra to epochs -20d < t < +20d, we restrict its application to SNe 2002gr and 2002go 
(we plot the solution corresponding to the smallest^ 2 - see text.), for which we have not plotted the gaussextract output for sake 
of clarity. The symbol © and arrows have the same meaning as in Fig. ^] Lower panel: Ratios of the various spectra to the 
specinholucy output. We have deliberately restricted the y-axis range to [0,2] to be able to visualise small ~ 20% differences in 
the SN flux. 



M(/t rest Xz,a,j8) = aS0Wl 



f/^OWl+z)) (a) 
z)) \ £S0W (b), 



of the supernova type and redshift is greatly improved upon and 
the en hanced accuracy of this approach is discussed in lSaintonl 



where A KSt and /l b s are the rest-frame and observed wave- 
length, respectively, (a) corresponds to the case where a local 
galaxy template spectrum is redshifted to fit the observed spec- 
trum, whereas (b) makes use of the observed SN host galaxy 
itself. The^ 2 is made robust against aberrant points (remaining 
sky lines, mostly). The fit is done for every galaxy-supernova 
pair and the results are sorted in ascending^ 2 . One is therefore 
able to evaluate how significant is the best result with respect to 
other SN-galaxy combinations, and in turn how confident one 
should be about the output supernova type, redshift and epoch. 

It should be noted that, contrary to the other methods pre- 
sented in this section, SN-fit is a postprocessing tool and thus 
cannot be required to correct for systematic errors made in the 
reduction process (most likely using standard software pack- 
ages such as IRAF or MIDAS 6 ), such as improper removal of 
sky lines. These errors will propagate in the SN-fit result, as 
seen in the output spectrum of SN 2002go where the Si II fea- 
ture is also affected by the atmospheric A-band. SN 2002gr on 
the other hand illustrates the limits of using a^ 2 minimisation 
technique to separate the SN from the underlying background 
in low signal-to-noise cases. 

The advantages of this purely statistical approach is that no 
a priori assumption is made about the underlying background. 
The host galaxy contamination is removed using a galaxy tem- 
plate or better still the SN host galaxy itself. The determination 

6 ESO-MIDAS is a registered trademark of the European Southern 
Observatory; see http://www.eso.org/projects/esomidas/ 



( 2004). The disadvantage is that we are limited by the database 
of local SNe, and any spectral peculiarities will not be inde- 
pendently accounted for, if they are not reproduced in local SN 
spectra. Also, the fact that in both cases we use the observed SN 
host galaxy itself as a template to decompose the extracted SN 
trace into a supernova and a galaxy component - case (b) in the 
above system of equations - illustrates what we were implying 
in section @] about the propagation of supernova "features" in 
the final SN spectrum due to the inherent imperfections of the 
galaxy template removal approach. 

8. Conclusion 

We have presented the advantages of using a more elaborate 
spectral extraction algorithm for obtaining supernova spectra 
devoid of background contamination, and hence enable a tru- 
ely quantitative analysis of high-z Type la supernova spectra. 
This will allow us in the future to set constraints on potential 
evolutionary effects that could affect SNe la as function of red- 
shift. Not only are we able to extract the PSF-like components 
of an input spectrum, one can also restore the entire 2D spec- 
trum and hence get a hold on systematic errors linked with the 
restoration process. We thus have a quantitative way of judging 
whether or not the extracted SN spectrum is optimal for sub- 
sequent analysis and whether we are not mistakenly "seeing" 
SN-like features in a pure galaxy spectrum, as is sometimes 
the case in low S/N high-z SN la spectra. 

However, there are some caveats mainly associated with the 
complexity of specinholucy. Errors associated with the use of 
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an inappropriate SSF for restoring the point-source channel or a 
mis-evaluation of the width of the spatial resolution kernel nec- 
essary for the point-source/extended background discrimina- 
tion will strongly affect the restoration (see section^}. Through 
restoration of the 2D background spectrum one is however able 
to produce diagnostic plots to critically evaluate the error in 
each case and repeat the restoration with different settings (see 
Fig.0. 

Note that a valuable by-product of this novel technique is 
the acquisition of a host galaxy spectrum with no contamina- 
tion from the supernova flux! Galaxies at z ~ 0.5 are at most 
a few arcseconds in size and one can essentially obtain inte- 
grated spectra of high-z SN host galaxies by centering a slit on 
both the SN and its host. This in turn enables the study of prop- 
erties of high-z SN la host galaxies and compare them with 
integrated spectra of local SN la hosts ^Sullivan et aTll2003t 
Willi ams et all2003tlGh allager & Garnavichl2004 . 

Thus, specinholucy provides the high-z supernova commu- 
nity with a reliable tool to obtain clean SN la spectra and in turn 
to derive accurate physical quantities associated with the ejecta. 
Two of the authors (Stephane Blondin and Bruno Leibundgut) 
are part of the ESSENCE project, and Gregory Sainton is part 
of the SNLS collaboration. Both teams recognise the need for 
a quantitative analysis of their results in order to constrain po- 
tential systematic effects that could distort the cosmological 
results derived from observations of high-z SNe la. We have 
shown specinholucy to be the ideal tool for this purpose. In fact, 
the VET data taken for the ESSENCE proj ect will be shown in 
their restored form via specinholucy in lMatheson et all2005[ 
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